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We study the equation of state (EOS) of nuclear matter as function of density. We expand the 
energy per particle (E/A) of symmetric infinite nuclear matter in powers of the density to take into 
account 2,3,. . .,N-body forces. New EOS are proposed by fitting ground state properties of nuclear 
matter (binding energy, compressibility and pressure) and assuming that at high densities a second 
order phase transition to the Quark Gluon Plasma (QGP) occurs. The latter phase transition is due 
to symmetry breaking at high density from nuclear matter (locally color white) to the QGP (globally 
color white). In the simplest implementation of a second order phase transition we calculate the 
critical exponent S by using Landau's theory of phase transition. We find <5 = 3. Refining the 
properties of the EOS near the critical point gives 8 — 5 in agreement with experimental results. 
We also discuss some scenarios for the EOS at finite temperatures. 
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I. INTRODUCTION 



In recent years, the availability of new heavy-ion accelerators which are capable of accelerating ions from a few 
MeV/nucleon to GeV/nucleon has fueled a new field of research loosely referred to as Nuclear Fragmentation. The 
characteristics of the fragments produced depend on the beam energy and the target-projectile combinations which 
can be externally controlled. Depending on the beam energy, hard photons, pions, kaons and so on can be produced as 
well [I |, ■ Fragmentation experiments could provide informations about the nuclear matter properties and constrain 
the EOS of nuclear matter. A 'conventional' EOS provides only limited information about the nuclear matter: the 
static thermal equilibrium properties^. In heavy ion collisions non-equilibrium processes are very important, thus 
nuclear transport properties will play an equally important role. If we want to study the EOS at high densities and 
high temperatures, we have to rely on theoretical estimates as well. The low density behavior of nuclear matter 
determines the observables mechanism of the final expansion stage in a collision before the break up. The low density 
part of the nuclear EOS is directly related to the final fragmentation, nuclear compressibility, momentum dependence, 
etc After an energetic nucleus-nucleus collision in the 100 MeV-4 GeV/nucleon beam energy region, many light 
nuclear fragments, a few heavy fragments and a few mesons (mainly pions) are observed. Thus the initial kinetic 
energy of the projectile leads to the destruction of the ground state of nuclear matter and converts it into dilute gas 
of fragments which then loses thermal contact during the break-up or freeze-out stage. 

One of the standard methods to explore the nuclear EOS is within the framework of the mean field theory. It starts 
out with a Langrangian including the nucleoli field a scalar meson field (f>, and a vector meson field V^. Customarily 
the contribution of the scalar field is described by a quartic polynomial. 

From conventional nuclear physics we know that there is a stable equilibrium state at the normal nuclear density 
po = 0.145 — 0.17/m _3 [3, H with a compressibility in the range of K = 180 — 240 MeV @ and a binding energy of 
15-16 MeV/nucleonQ. In this work we take po = 0.165/ra 3 and K = 225 MeV based on the condition that the mean 
field potential has a minimum at normal nuclear density. With increasing density, the effects of N-body correlations 
become more and more important. This is especially true near a phase transition. Furthermore, nucleons are not 
elementary particles but they are made of quarks and gluons, thus N-body forces are expected to be stronger at high 
densities where the nucleon wave functions strongly overlap. We can take these features into account by expanding 
the EOS in powers of the density as customary in the virial expansion of any EOSQ. 

Finally we discuss some properties at finite temperatures assuming either a classical gas or a quantum Fermi 
systems. We show that at the densities and temperatures of interest the classical approximation is not valid. This is 
at variance with many experimental and theoretical results in heavy ion collisions near the Fermi energy 0-[ll| which 
assume the classical approximation to be valid. Quantum corrections have been recently extensively discussed in (l2| 
for cluster formation in a low density expanding nuclear system. 



II. NUCLEAR EQUATION OF STATE 



For a system interacting through two body forces having a short-range repulsion and a longer-range attraction the 
EOS resembles a Van Der Waals one. This is indeed the case for nuclear matter [ll43l. Il3j. 

A popular approach is to postulate an equation of state which satisfies known properties of nuclei [7| . The equation 
for energy per particle is: 



E/A = 22.5pi + -~p + — -F (1) 

Z 17+1 

where p = p/ po and po is the normal nuclear density. The first term of eq. (TTJ) refers to the kinetic energy of a Fermi 
gas. The other terms are due to potential interactions and correlations. To generalize at finite temperature we could 
use a classical approximation giving an EOS in the form [l3j . 

P = „« + „r (2, 

In the following we will test the validity of such an approximation. 

The requirement of causality provides several theoretical constraints on the EOS [14] at high densities and limits 
the choice of the functional form of the compressional energy that can be used inphenomenological EOS (see Fig. [4}. 
Very stiff equations of state may lead to superluminal speed of sound (see ref. [15|]). However, this is not a problem 
if the acausality occurs in a region of the phase diagram where matter is in the mixed or plasma phase, because the 
phase transition softens the EOS. There are some phenomenological parameterizations of the specific energies like 
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Linear Quadratic, Sierk-Nix or Grant-Kapusta which are acausal at sufficiently high densities Fortunately the 
acausality occurs well within or beyond the mixed or plasma phase for all the parameterization except the Quadratic. 
The basic parameter is the (isothermal) compressibility, which is defined as 



dP 
dp 



(3) 

p=p o ,T=0 



The three parameters A, B and a in eq. {T]) are determinated by using the conditions of pressure equal zero, 
the binding energy of E/A = — 15MeV and finally the compressibility is of order of 200MeV (as inferred from the 
vibrational frequency of the giant monopole resonance[6]), at p — pq. Using these conditions, we get A = — 356MeV, 
B = 303MeV and a = 7/6, we will refer to this EOS as CK200 (Conventional, K=200 MeV). 

Now if we modify this approach in accordance to [16] the compressibility condition is substituted by the mean held 
potential having a minimum at the ground state density: 

U(p) = A(p) + BA{~pY (4) 

where 

ou(p) 



dp 



(5) 



P=Po 

which means that the compressibility at p = po is equal to K = 225 MeV then we have the conditions 

a) E/A = -15MeV 

b) P = at p = p (6) 

c) K = 225MeV 

Solving these equations, we get A = -210 MeV, B = 157.5 MeV and a = 4/3, we will refer to this EOS as CK225. 

The form of the EOS is a very delicate subject. For a nuclear system we expect to see a liquid-gas (LG) phase 
transition at a temperature of the order of 10 MeV and at low density. Under these conditions we can assume that 
nuclear matter behaves like a classical ideal gas eq. ([2]) however this is just our ansatz. Actually in this work we will 
compare the theoretical behavior of the EOS assuming classical ideal gas and Fermi gas. 

In order to calculate the critical point, we will impose the conditions that the first and second derivative of eq. @ 
respect to density are equal to zero therefore we can obtain the critical temperature and density, i.e. T c = 9 MeV and 
p c = 0.3/?o- Such values are in some agreement with experimental results (l_3j|. However, we notice that in order for 
the classical approximation to be valid, the ratio of the temperature to the Fermi Energy, e/(p) = 36(p/po) 2 / 3 MeV, 
should be much larger than one. For the values above we get T c /ef — 9/16 < 1 which shows that we are still in a 
quantum regime. 

The validity of eq.(l) is restricted to densities close to the ground state value. In fact no further constraints are 
imposed so far for larger densities. Such constraints should come from experimental data in heavy ion collisions and 
properties of heavy stars. Those data, if available, do not give directly a constraint on the EOS but must be filtered 
through model calculations. The models in turn need some form of EOS. We propose a new equation for the energy 
per particle which could be used in microscopic calculations: 

k 

E/A = e f pi + £ (7) 

n— 1 

where the first term refers to the average kinetic energy of a free Fermi gas with if = 3/5e/ = 22.5MeV, the 
other terms are due to potential interactions and correlations. The n = 1 term is obtained by taking into account 
the interaction between pair of particles, and the subsequent terms must involve the interactions between groups of 
three, four, etc., particles. The coefficients A n in the expansion, eq. ([7]) are called first, second, third, etc., virial 
coefficients @. 

Let's start considering three body forces 0(p 3 ) and we assume that at the normal density and zero temperature the 
energy of the ground state is -15 MeV and the pressure is zero. In addition we assume that the mean field potential 
has a minimum at normal density or, equivalently K — 225 MeV. If we do this, we have three conditions and three 
equations so we can solve the corresponding set of equations. Unfortunately, the solution has no physical meaning 
because the energy diverges to minus infinity when the density approaches infinity, see table I. 

For the fourth order of our expansion the eq. ([7]) takes the form 

~ , . 2 t4i A<2 o A% q Aa a . . 

Ej A = ejp* + ^P+ ~fp 2 + ~fp 3 + -fp^ (8) 
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Figure 1. The energy per particle of nuclear matter as a function of density using different density-dependent interactions in 
comparison with the 'conventional' formulation. The parameter values are given in Table [I] Symbols: CK225-thin solid line, 
CK200-dashed line, CCS<53-dashed dotted line and CCS<55-thick line. 



Table I. Values of the parameters of the EOS's of nuclear matter. 



Interactions 


A 1 

[McV] 


A 2 

[McV] 


A 3 

[McV] 


A 4 

[McV] 


A 5 

[McV] 


A a 

[McV] 


a 


CK200 


-356 


303 










7/6 


CK225 


-210 


157.5 










4/3 


3 


-135 


112.5 


-30 










CCS53 


-136.89 


120.99 


-41.32 


4.72 








5 


-137.59 


124.41 


-46.51 


7.67 


-0.47 






CCS55 


-137.96 


126.25 


-49.46 


9.55 


-0.92 


0.0035 





The 'Conventional' EOS corresponds to eq. Q where the the ground 
state binding energy E m i n =-15 MeV and the compressibility at normal 
nuclear density K=200-225 MeV and the remaining EOS's correspond to 
eq. (0 with E min =-15 MeV and K=225 MeV respectively. 



Here we assume the symmetry breaking at high density, from nuclear matter (locally color white) to the QGP 
(globally color white), this possibly gives a second-order phase transition (l7l [Tc| . 

In accordance to the conditions given in eq. ^ we can add two extra constraints based on the conditions of matter 
close to the critical density of a second-order phase transition at T=0. 

At the critical point the first and second derivative of the pressure with respect to density are equal to zero. Then 
we have five constrains: 
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a) E/A = -15MeV 

b) P = at 

c) K = 225McV 

«0 f | =° 



P = Pa 



(9) 



Solving these equations we get the values of A\..A^ (see table |TJ) . Therefore we obtain the critical point p c — 2.9354, 
we refer to this EOS as CCS<53. This could be the critical point for a second-order phase transition to the QGP at 
T=0. 



In order to include higher order terms we assume that 



= where n — 1, ... ,4 (this assumption will 



P=Pa 



become clear later on) then we solve the resulting nonlinear system of many unknown variables to get all A n and 
therefore the critical point p c at T=0. Here we are taking into account the interaction between pairs of particles and 
interactions between groups of three, four, five and six particles. 



For interactions up to five particles 0(p 5 ) where 8 ' 



dp" 



with 11=1-3, the nuclear EOS diverges negatively 



P=Pc 



again as in the 0(p 3 ) case. 

If we take more variables in our expansion, up to sixth order we get: 



Ej A = iff)* 



A 



A 



3 ~ 3 



A 



V 



A 



6 ~6 



(10) 



Imposing the conditions that the fourth order derivative of the pressure vanishes as well, gives the values of the 
parameters reported in Table I and a critical density p c = 5.2, which we refer as CCS<55. 

Now we can try to apply our equations of state in the high-density domain. A quite promising result is shown in 
Fig. [TJ where we compare the different EOS at T=0. The different EOS are very similar near and below the ground 
state density of nuclear matter while they differ greatly, as expected at higher densities. In particular the EOS softens 
because of the assumed QGP phase transition. 

The pressure is shown in Fig. [5] at temperature T=0 with a comparison with the conventional EOS. At p = po the 
pressure is zero and increases largely for the 'conventional' EOS. 

We proceed now to the study of the compressibility. We are assuming that at p = po, K=225 MeV (apart CK200) 
as is shown in Fig. [3] A negative compressibility (which indicates an instability region) is obtained at subnuclear 
densities only, while it becomes zero at the critical point for the QGP phase transition. 

A negative compressibility gives an imaginary speed of sound since for a particle of mass m [19j: 




(11) 

Such a quantity is plotted in Fig. 4 as function of density for the different EOS. Two phenomena are worth noticing. 
First the speed of sound becomes larger than the speed of light for all the EOS (excluding CCS^5) for p > 4.5po- The 
CCS<55 EOS gives a superluminal speed of sound at almost twice such a density, well in the region of the QGP. 

The second property is that the speed of sound becomes imaginary in the instability region, thus a discontinuity 
is shown in Fig. 4 below normal nuclear matter density. The small discontinuity observed near the critical density 
of CCS<55 is due to the numerical solution of the set of equations used to determine the values of the coefficients 
reported in Table I. We can have a better view of the critical region by zooming on the density as in Fig. 5. Now 
the discontinuities at lower densities due to the LG phase transition are visible for three EOS. The instability region 
is very much the same for all the EOS having the same compressibility, which shows that such a region is mainly 
determined by the ground state properties of the EOS and not by the assumed functional form. As we will show later 
also the critical point for the LG is the same when the compressibility is the same. At higher densities, the speed of 
sound becomes zero (discontinuous) at the critical densities of the fourth and sixth order EOS respectively. 
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Figure 2. The pressure per particle of nuclear matter as a function of density using different density-dependent interactions. 
Symbols: CK225-thin line, CCS53-dashed dotted line and CCS<55-thick line. 

III. FINITE TEMPERATURES 

In order to study the properties of the EOS at Unite temperatures we need to go b eyond the classical approximation. 
A simple functional form could be obtained using a Fermi gas expansion^, HH Il3j . 

The region of validity of the Fermi gas model is related to the ratio of the temperature to the Fermi energy T/ef(p). 
In this scenario eq. ([7]) takes the form: 

k A 

Sin 



E/A = i f pi + J2 -f?P" + *oP^T 2 (12) 
where the level density parameter clq = 1/13.3 MeV -1 . Correspondingly the pressure takes the formQ: 



P = Po 



5 A ni n p( n+1 ) 2 i , 
^ 3 



(13) 



The properties near the critical temperature and density can be obtained for the various EOS imposing the con- 
straints similar to those given in eq.(9d),e)). For the case of fourth order 0(p A ) the critical temperature and critical 
density for the LG phase transition has the values T c = 18.05 MeV and p c = 0.3724po respectively. In this approach 
we are assuming interactions between pairs of particles and between groups of three and four particles which could 
be fully justified by the fact that nucleons are made of quarks and gluons. For higher order terms, 0(p e ) the critical 
temperature and density are T c — 17.932 MeV and p c = 0.3717po- Similar values are obtained for the CK225 EOS. We 
notice that the critical values of the LG do not change much by using different EOS. In particular the critical density 
seems to be almost independent of the assumed classical or quantum statistical properties of nuclear matter. On the 
other hand, the critical T changes almost of a factor two when going from the classical to the quantum approximation. 
We have already noticed that the ratio of T to the corresponding Fermi energy at the critical density was smaller 
than one when using the classical approximation. If we perform the same calculation for the quantum case we get 
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Figure 3. The compressibility per particle of nuclear matter as a function of density using different density-dependent interac- 
tions. Symbols as in Fig. 2. 



T c /ef — 18/18.6 ~ 1, which implies that higher order terms must be considered in the expansion of the Fermi gas 
pressure (or energy per particle) but we are not in the classical regime. A detailed calculation at finite T properly 
taking into account the Fermi statistics is beyond the scope of this work but we would expect a slight reduction 
of the temperature as compared to the value given above. We notice again that values obtained experimentally for 
the critical T and p sug gest that we are in the quantum regime while the methods used to obtain them are purely 
classical [9l-flll|. but |l2| This is analogous to what we have done above, eq.(2), using a classical approximation for 
finite T. We expect that properly taking into account quantum statistics will dramatically change the values obtained 
experimentally near the critical point. 

In the quantum approximation considered above we can obtain the P(p,T), the result is given in Fig. 6 for the 
CCS<55 case. This region is very similar for all the EOS for a given compressibility (K=225 MeV in our case). At 
higher T and/or p the pressure flattens because of the QGP phase transition and increases again at higher densities. 

We can compare the differences between different EOS more in detail in the purely classical, Fig[7] and quantum 
approximation, Fig[8l We notice the dramatic differences between the two cases especially at the lowest densities and 
T where quantum effects are stronger. For completeness we show in Figs. 9 and 10 the same cases but as function of 
reduced volume. 



IV. CRITICAL PHENOMENA 



The term critical phenomena refers to the thermodynamic behavior of a system near the critical temperature of 
a second order phase transition. A simple understanding of this phenomenon can be obtained in the framework of 
Landau's theory [8]. In such a framework it is assumed that near the critical point of a second order phase transition, 
the relevant degree(s) of freedom reduce to a few (order parameters) which reflect basic invariance properties of 
the Hamiltonian. Below the critical point such invariance(s) is spontaneously broken and it is restored above the 
critical point, this means that the order parameter is non zero only below the critical point. Some examples are the 
magnetization in a ferromagnetic system, or the difference in densities between a gas and its liquid in a liquid-gas 
phase transition. In Landau's approach one assumes that there exist a free energy 'F near the critical point that 



9 




40 



Figure 6. Equation of state surface for a nuclear system with a second-order phase transition to the QGP, CCS55. The pressure 
is reduced by the phase transition near the critical region. The pressure increases again at very high densities in the QGP 
phase. 



depends on the order parameter and its conjugate field. In general "J" = E/A — TS, where S is the entropy. Near the 
critical point we can expand the free energy in terms of the density at the critical temperature (which is zero for our 
case). Defining the order parameter r/ = V — V c i.e. the distance from the critical volume, we get 

(V-V c )^ (14) 

v c 

n=0,l...N. We define an external or conjugate field P(V C ) = —dE /dV\v=v c - If we stop our expansion at fourth order 
n = 4. and notice that in the neighborhood of the critical point the conditions of the second order phase transition are 
P'(V C ) = and P"(V C ) = at T = . Thus only n=0 (a constant), 1 and 4 survive in the expansion above. Imposing 
the free energy to have a minimum in presence of a conjugate field gives 



rj ~ pi ~ pi (15) 

Therefore 0(p 4 ) has the critical exponent 5 — 3 which corresponds to the 'mean field' or 'classical' valueQ. Exper- 
imentally the value of S = 4 — 5[8( is found. It is now evident how to go beyond the 'mean field' value using the 
arguments given above. In particular if we impose P"'(V C ) = and continue the expansion in eq.(14) to n = 5 we 
easily get 8 = 4. However the resulting fifth order coefficient in the virial expansion of the EOS is negative, see table 
I, thus this case is unphysical. Finally expanding to n = 6 and imposing P""(V C ) = gives 8 = 5 and the fitting 
parameters reported in table I for the CCS(55 EOS. Notice again that our formulation is perfectly consistent with 
Landau's theory as discussed for general thermodynamical systems. Because of this analogy we expect that all other 
critical exponents are the same as those calculated in the Landau's approach [8j. 



_ 1 d {n) E(V) 
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Figure 7. The pressure per particle of nuclear matter as a function of density at some temperatures. Here we are assuming that 
nuclear matter behaves like a classical ideal gas. The EOS are respectively from top to bottom: CK225, CCS<53 and CCS<55. 

V. SUMMARY AND CONCLUSIONS 

In this paper we discussed some nuclear equations of state bridging basic properties of nuclear matter near its 
ground state and a possible second order phase transition to the quark-gluon plasma at high densities. We have 
determined a critical density of about three times normal ground state density for the QGP in the case where a 
critical value S = 3, this corresponds to the so called mean field value or classical value of the critical exponentQ. 
To go beyond this classical value we used Landau's theory of phase transition and for <5 = 5 we determined an 
EOS to 0(p 6 ) which gives a critical density about five times the ground state density of nuclear matter. Using an 
MIT bag model, it is possible to estimate a critical density for the QGP of about five times normal nuclear matter 
density using a bag constant B 1 / 4 = 206MeV[I3. Other estimates using phenomenological hadronic and QGP EOS 
give similar densities but a first order phase transition[3]. In order to see if there is a first or second order phase 
transition or a simple cross-over from one state to the other, reliable experimental data in the region at high baryonic 
densities and relatively small temperatures are needed or data from heavy stars compared to refined theoretical 
models. Such quantities could be for instance collective flow compared to microscopic calculations which implement 
the EOS discussed here. As we have seen the pressure at high densities is completely different in the EOS with or 
without a QGP phase transition. In particular if the transition is second order than the corresponding EOS cannot be 
much different from the one estimated here when S — 5. At lower densities a liquid gas phase transition might occur 
at the critical density p c = 0.37/?o and T c — 18 MeV. Such values are rather independent from the EOS chosen when 
fixing the compressibility. However the critical temperature depends somewhat on the assumed classical or quantum 
statistics. In the region of densities and temperature of the LG transition it seems that there is no reason why a 
classical approximation should work. However, presently most experimental and theoretical results in this region are 
obtained using classical methods. This is a feature that should be improved in future works in order to have more 
reliable constraints on the EOS|12|. 

Finally, we found that for odd orders as, 0(p 3 ), 0(p 5 ) this approach is not suitable to describe the basic properties of 
our system. The even order 0(p i ), 0(p 6 ) approximation result in more acceptable behavior of the energy per particle 
as function of density. It could be possible that the collective character in regions with high density of particles are 
associated with even numbers [2(1 . 
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Figure 8. Same as Fig. 7 but for a Fermi gas. 
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Figure 9. Behavior near to the critical point for 0(p 4 ); the pressure per particle of nuclear matter as a function of volume at 
some temperatures. Here we are assuming that nuclear matter behaves like a Fermi gas. 
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Figure 10. Behavior near to the critical point for 0(p e ); the pressure per particle of nuclear matter as a function of volume at 
some temperatures. Here we are assuming that nuclear matter behaves like a Fermi gas. 
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